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Abstract 



The problem in identifying the sites of origin of Galactic Cosmic Rays (CRs) is reviewed. 
Recent observational evidence from very-high energy (VHE, energies above 100 GcV) 
7-ray measurements is in contradiction with the surmise that synchrotron radiation 
from relativistic electrons is indicative for hadron acceleration. It rather points to a 



less than the value required if these objects were to be the main sources of Galactic 
cosmic rays (about 30%). Observations of CR anisotropy and the emission of low- 
energy (energy < 10 GeV) 7-rays from the Galactic disk indicate that the sources of 
low-energy cosmic rays are distributed with a Galactocentric radial scale length in the 
order of 25 ± 10 kpc, much larger than expected if SNRs are the main sources of CRs. 
These two facts - together with the body of evidence from CR isotope abundances - 
strongly suggest that a new class of astrophysical objects - distinct from SNRs and 
located manly in the outer reaches of our Galaxy - is the major source of hadronic CRs 
in our Galaxy. 

The basic observational features of ultra high-energy (UHE, energy > 10 19 eV) CRs 
are most naturally understood if the same CR sources accelerate CRs up to the highest 
observed CR energies. Proposals for the nature of a new source class are mentioned. 
The origin of CRs is still as much shrouded in mystery as it was in 1957, when Philip 
Morrison wrote a seminal review about CR origin. The potential for discoveries is thus 
great. 



CR-acceleration efficiency of supernova remnants 




below one percent, much 
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1 The term "SNR" is below always used in the sense: "the well known non-relativistic ejecta of supernova 
with a typical kinetic energy of 10 51 ergs" 



Thus in the 40 plus years since the publication of Morrison's seminal paper^fq] ... we have 
come no closer to a definitive model of cosmic-rays origin. 

Trevor Wockos, August 1999 |lisj ] ■ 

Until recently the hypothesis that cosmic rays originate in the flashes of supernovae has 
been based solely on energy considerations - evidence that is still far from adequate. 

Vitaly Ginzburg, July 1953[Ba] • 



1 Introduction 

1.1 The problem in historical perspective - outline 



Around 1930 Millikan,Cameron[73] and Regener pql realized that the energy density of cos- 
mic rays in space is about as high as that of integrated star-light. This made clear that 
cosmic rays are an important cosmic phenomenon, and two questions were raised. 

1. By which mechanism are cosmic rays accelerated? 

2. Where and in which sources does this occur? 

In this manuscript I try to summarize what is presently known about the answer to the 
second question. I strive for brevity and will review neither the theoretical nor the ex- 
perimental literature exhaustively. My style is inspired by a memorable article of Philip 



Morrison with almost the same title in 1957 [76] ("On the Origins of Cosmic Rays" ), be- 



cause our situation concerning the second question is hardly better than the one described 
by Morrison like this: "the broad road (towards the answer) is finally visible, but the details 



may well be all wrong\76}" . 

Up to 1948 both questions remained a complete riddle, revealing a deep fundamental lack 
of knowledge of the physical world. Highly speculative origin theories - without a firm 
basis in contemporary physics - were proposed. Millikan|74T| suggested a new fundamental 
production process and Lemaitre p^j explained cosmic rays as relics of the early universe. 
These ideas stood at the cradle of modern nuclear astrophysics and cosmology but did not 
really advance cosmic-ray physics. Two publications by Baade and Zwicky from 1934|J, 
suggesting an origin of CRs in extragalactical supernovae, stand out for their impressing 
far-sight|. 

Then, in a brief, productive "golden period" - from about 1948 to 1962 - Fermi conceived 
the fundamental theoretical concept for the acceleration mechanism and several crucial ex- 
perimental discoveries provided the basic framework in which we still grapple today for an 
answer to the question of the site of origin. Fermi's brilliantly simple idea was the real- 
ization that nuclei colliding via magnetic interactions with macroscopic gas clouds reach 

2 Their basic idea was that relativistic particles with a total energy between 10 53 and 10 54 ergs would be 
ejected in a supernova. These particles indirectly "heat" the stellar debris and thus lead to the visible SN 
event. It might still turn out that this image is very close to being correct. 
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equilibrium only when their energies are similar, i.e. when the nuclei possess macroscopic 
kinetic energies. Henceforth the question was no longer "how can such high energies be 
reached?" but "under which circumstances and with what time scale?" This was the death- 
blow for the early speculative theories, CR origin had become a problem of conventional 
astrophysics. Since the mid 1980s very speculative theories are considered once again, see 



section 4.2 



The most important experimental achievement of the "golden period" was arguably the 
determination of the CR spectrum as a function of total particle energy. Today the flux 
of the dominating hadronic component, consisting of protons and heavier nuclei, of the 
local non-solar cosmic rays (CRs), is determined at energies between about 0.1 GeV and 
3 • 10 20 eV. Neglecting solar modulation effects, the all-particle CR flux can be described 
by power-laws with indices in the narrow range between -2.5 and -3.2 over the whole range 
(more than 10 decades) in energy (see data points in fig. [l]) and a single power law with an 
index of -3 is quite a good approximation to the experimental data (see thin dashed line in 

A framework for the explanation of this spectrum, assuming various different sites of ori- 
gin, is at present practically universally accepted and is sketched in fig]l[ In sections |2| 
and ^| I will critically discuss the origin of the low-energy and very- high energy component, 
respectively. Following a seminal paper by Biermann and Davis section ^ explains why ob- 
servations make necessary a departure from the conventional framework. Section |B| briefly 
discusses two recent and one historic proposed alternative major source of CRs. In the final 
section |6] I will give a personal interpretation of the current situation and lay down some 
questions to be answered by future research. 

1.2 Cosmic-Ray Origin: the single most important problem of "As- 
troparticle Physics" 

When I talk to younger colleagues in the field of astroparticle physics I get sometimes the 
impression that for them, the question of low-energy cosmic-ray origin is some old problem 
left over from the last century. Nothing could be further from the truth. Cosmic-ray origin 
is the only important physics problem where astroparticle physics is called upon to take the 
lead - both theoretically and experimentally - in providing a full solution for all energies. 
In other areas, like the physics of active-galactic nuclei, the emission-mechanism of pulsars, 
the quest for dark matter or the nature of cosmological background fields astroparticle 
physics just supplies some pieces to jigsaw puzzles being assembled by others. 

2 Origin of "low-energy" (energies below the "knee") CRs 
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2.1 The proposal of SNRs as sources of the low-energy component 

In 1953 Iosif Shklovsky ||9~5| - armed with his recent identification of the bright radio source 
Cas A as the remnant of a historic supernova explosion (SNR) f^] - proposed that SNRs 
(rather than SNe themselves) are the main sources of Galactic cosmic raysQ. Shklovsky 
gave two arguments for his idea: 

1. the total power of the of CRs injected into the Galaxy is "similar" to the kinetic explo- 
sion power of a Galactic supernovae (i.e. 10 51 erg/50 years) []. Quantitatively the power of 
the injected CRs was determined to be about 30 % (with an uncertainty of a factor 3)[f47|] 
[] of the total kinetic power of young Galactic SNRs[]. 

2. SNRs are most prominent non-thermal radio sources because they produce non-thermal 
electrons copiously. Shklovsky wrote |)5||: "...it is quite natural to consider that simultane- 
ously with relativistic electrons, also occur relativistic heavy particles^ 



2.2 A guiding principle: is synchrotron emission indicative of the total 
non-thermal acceleration power? 

Shklovsky's second point fell on fertile ground. In the year before Ludwig Biermann had 
made the "tentative assumption that the radio emission indicates the order of magnitude of 



all non thermal emissions. "|16j. Henceforth, the assumption that the main CR sources are 
prominent radio sources was a beacon to most theoretical studies on CR origin. In 1962 
Geoffrey Burbidge formulated the idea that the other class of celestial radio sources that 
rival SNRs in apparent luminosity, active galaxies (typified by Cyg A, the second brightest 
radio source in the sky after Cas A) are the main sources of all cosmic rays J 



The following decade's discussion concentrated on a heated controversy between Burbidge 



and Vitaly Ginzburg, who had become the main proponent of Shklovsky's idea|54, 55]. A 
discussion of alternatives (such as e.g. a stellar origin of a very-low energy component as 
favored by Morrison [76] or the ideas of Biermann and Davis [17], discussed in section |3|) fell 



completely by the wayside. 

3 An identification that was initially vehemently resisted by Walter Baade[Q] who had proposed SNe as 
cosmic-ray sources some 20 years earlier. 

4 Shklovsky's priority on this crucial idea is clearly stated by V. Ginzburg in an early review [|53f . 

5 This "equality" does not point directly to SNRs as the site of CR acceleration, rather it only points 
towards a possible connection "CRs - supernovae". 

This is the most recent value (a previous careful determination gave the range of 10 - 30 %Q). In 
Shklovsky's article the total injected power by SNRs was estimated to be about 2 orders of magnitude 
smaller (this was sufficient given the incomplete knowledge of Galactic CR propagation in 1953). Such a 



smaller efficiency is in line with the expectation of Biermann and Davis (see section 3.1). 

7 Such an argument was first made made by Baade and ZwickyQ in an extragalactic context. With the 
same values and assumptions about propagation and explosion energy as Shklovsky, it was first made by 
ter Haar[|io4|. 
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2.3 The present situation: are SNRs the main source of low energy CRs? 

In spite of enormous advances in experimental techniques in the half century that has 
elapsed since the publication of Shklovsky's paper, not a single additional observa- 
tional fact directly supporting the hypothesis that SNRs are the major source 
of Galactic cosmic rays has come to light. 

To make matters worse - in what I consider to be the most important experimental discov- 
ery in CR physics since the 1950s - it was recently learnt that Cas A is not an important 
source of hadronic cosmic rays (see section |3.1| for further discussion) . This saps the foun- 
dations of Shklovsky's second argument. The beacon mentioned in the previous subsection 
- that still guides many0 - threatens to have misled us for nearly half a century. The 
dominant sources of hadronic CR are either SNRs or they are not prominent 
radio sources. The sole direct evidence in favor of a SNR origin of CRs is again the ter 



Haar's |104f| energy argument, evidence that is still far from adequate" p2 \ 
But if not SNRs, then what^^ ? Let us take a big step back in time and reconsider some 
early arguments against SNRs as CR sources in the light of today's body of experimental 
evidence, in an attempt to point a way out of this impasse. 



3 A Path not taken - 1958: Returning to a critique by Bier- 
mann and Davis 

Near the end of the "golden period", in 1958, two of its principal architects, Ludwig Bier- 
mann and Leverett Davis, summarized some lessons learnt [|T^| . Their analysis, though 
certainly noticedP], had little impact on the consecutive research)^. This is a pity since their 
doubts and some of their surmises have been born out by observations in a truly impressive 
way. 



3.1 Acceleration efficiency in SNRs 

Biermann and Davis write: "With the present frequency of supernovae, it would just be 
possible to account for the acceleration of cosmic rays if one assumes a very high efficiency 
for the conversion of kinetic to cosmic ray energy. Even with a higher frequency in the 
early stages [of Galactic evolution], the efficiency remains so high (several percent) that the 
present authors find this proposal difficult to accept. " 



Today's value for the required efficiency at about 30 %[47] is even higher. Biermann and 



Davis refrain from telling us why to them an efficiency of several % seemed unrealistically 

8 The gist of a prominent colleague's only reaction to a recent manuscript |^| was : "Nobody believes 
that because the proposed CR-producing ejecta were not seen in the radio." 

9 Their collaboration on Galactic CRs was explicitely mentioned in the laudation on occasion of the award 
of the Bruce medal to Biermann jlj. 

10 The paper was quoted only five times, the last time in 1969. 
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high, but it is easy to guess the reason for this judgment. As pioneers of solar-system CR 
physics, they knew that the efficiency of hadron acceleration to relativistic energies in large 
solar flare events is typically 0.1 % fl6|f n | . The velocities, matter densities and magnetic field 
strengths involved in these events are not radically different to those in SN ejecta. However, 
there are also important differences in temporal and spatial scaldT^j so Biermann and David 
prudently did not cite solar-system evidence in support of their doubt. Today's direct evi- 
dence on the hadron acceleration efficiency in SNRs strongly supports their doubts. 
Accelerated hadrons interact with ambient matter producing ttq mesons. The level of 
very-high energy 7-ray radiation expected from ttq decay allows to derive the efficiency of 
proton acceleration! 42] and has been searched for in a considerable number of SNRs using 
air-shower arrays and air Cherenkov telescopes. Let's take a closer look at the two most 
interesting cases. 



3.1.1 SNR G78. 2+2.1 - a cosmic beam dump 

G78. 2+2.1 - probably the remnant of a core-collapse supernova in the Cyg OB9 association 
1.4 ± 0.6 x 10 4 years ago - is fourth brightest SNR at 1 GHz in the skyQ. At a distance 
of only 1.5 ± 0.3 kpc it subtends an angle of 1 degree, about twice as much as the moon. 
Based on various observational indications Pollock^] suggested that a bright "hot spot" 
(called "DR4") of radio synchrotron emission in the south-east corner of the SNRs radio 



shell is due to the dense molecular cloud Cong 8 [33], that was recently hit by the advancing 
SN shock front. Since the time Pollock wrote the paper, Fukui and Tatematsu[|50|] have 
reconfirmed a CO emission at 2.6 mm and Wendker et al. |114| found H2CO absorption in 
four out of five positions in DR4. Because CO and especially H2CO are indicators of high 
density, "...there is good evidence that interaction with quite a dense cloud is occurring in 
this part of the SNR. (cited from | {114fl )". The reality of the association between shock 
wave and molecular cloud has been doubted because Frail et al.|4S] failed to find OH(1720) 



MHz Maser emission from SNR G78+2.1. This is unwarranted, Frail et al.|4S] write:" .Since 
masers are a nonlinear physical phenomenon, depending sensitively on input parameters the 
non detection of OH (1720 MHz) maser emission does not in itself rule out an interaction." 
In fig.|| I compare the detailed shape of the IR emitting region at 100 fim as measured by 
the IRAS satellite (taken from IRAS ISSA server |6l[] ) and DR4 at 92 cm as measured by 
the Westerbork radio telescope (from the WENSS survey [|115| ). The position and general 
elliptical shape, with a major axis from SW to NE, of the IR and radio emitting region 

11 The stress is on relativistic, the efficiency for acceleration to lower, but still non-thermal energies can 
be much higher |Q. 

12 The acceleration efficiency for relativistic electrons in nova explosions - intermediate in energy between 
solar flares and supernova outbursts - has been been determined to about 1 %|9Cl]. This value is similar to 



the one in SNRs (section 3.3). Perhaps an efficiency for acceleration to relativistic energies in the order of 



0.1 - 1 % is typical for non-relativistic shocks in a magnetized low-density medium for all species. 
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are similar. This morphological resemblance suggests that indeed a molecular cloud was 
heated and compressed by the SNR shock wave. 

G78. 2+2.1 accelerates electrons to a spectrum with a power- law index of -2.08 ± 0.04, very 
close to the theoretically preferred value of -2,1 |116| ]. This makes G78. 2+2.1 an excellent 
candidate for a typical hadronic CR accelerator. The spectral index within DR4 is similar 



to the one in the rest of the remnant [116]. The entire cloud must then act as a "beam 
dump" for accelerated protons. The high matter density in the molecular cloud and rela- 
tive closeness of the remnant make DR4 a most sensitive test-bed of hadron acceleration in 
SNRs§. 

Fig.|] compares various observational upper limits on VHE emission from DR4 with theoret- 
ical expectation, calculated under the assumption that 10% of the kinetic explosion energy 
was converted to hadronic CRs. The most restrictive limit (from 47 hours of on-source 
measurement with the HEGRA system of telescopes [57]) falls below the expected value by 
more than a factor of 100. 



3.1.2 Cas A - cradle and bier of the belief that SNRs are the main source of 
CRs 

Cas A, at a distance of 2.8 - 3.5 kpc, is the remnant of Flamsteed's type lb supernova "3 
Cas" in the year 1680. It is the brightest non-thermal radio source in the sky and one 
of the best studied and understood non-thermal objects pO]. It was Cas A that inspired 



Shklovsky to propose that Galactic CRs are mainly accelerated in SNRs [ 95 1 . His idea was 
refined in the 1970s by models in which the CRs are mainly accelerated at the outer plane 
shock front via the first-order Fermi mechanism JTT], f41~| . Recently very clear experimental 
evidence from X-ray measurements for such a shock was found in Cas A (see fig. ||). 

Berezhko, Piihlhofer and Vork|i"3|] carefully calculated the expected VHE 7 flux from 
Cas A, under the assumption that CR acceleration efficiency of Cas A is 20 %. They 
point out that the observed flux of 7-rays with energies above 500 GeVQ falls below the 
expected flux by more than a factor of 100. The observed signal might well be due to inverse 
Compton emission by CR electrons, that are known to exist in Cas AQ. The discrepancy 
would then increase. 



3.1.3 SNRs are not the main sources of cosmic rays - their efficiency for ac- 
celeration is too low 

The upper limits on 7 radiation from these two remnants bound the hadronic accelera- 
tion efficiency to about 20/100=0.2 %. This shows that the intuition of Biermann/Davis 
was correct: the efficiency with which the kinetic energy of SNRs is converted to CRs is 
not much larger than 1 %, thus excluding SNRs as the major source of hadronic 
Galactic cosmic rays. 

Let us consider various attempts [0, O to avoid this conclusion. 
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The different ages of Cas A and G78. 2+2.1 exclude that the small VHE 7-ray flux is due to 
a long time scale of hadron acceleration ("late injection") or faster than usually expected 
diffusive CR loss from the remnant. 

There is a handful of other SNRs for which less restrictive or more model-dependent 



upper limits from Cherenkov telescope |10q, 22, t)2| and low energy-threshold air-shower 



arrays|75, 84] below the theoretical predictions exist. Therefore speculating that only a 
small fraction of "special" SNRs are the main sources of CRs quickly runs afoul with the 
required efficiency of particle acceleration that can hardly be much higher than the 30 % 
required if all Galactic SNRs contribute equally^ 3 ]. 

The non-detection of 7-rays in the required amount might indeed be due to an upper energy 
cutoff below the threshold of the detectors (at 200 GeV and higher for various Cherenkov 



telescopes and at 3 TeV and higher for various arrays). Berezhko and V61k[14] considered 
this possibility and aptly concluded: "Note however, that in this case SNRs would hardly 
be considered as the sources of Galactic CRs. Pf 7 

Finally, the enormity of the discrepancy for the two cases presented above precludes to 
invoke a SNR source spectrum steeper than the one expected from the standard theory 
(power-law index of -2.1). The steepest index compatible (with difficulty |70t| ) with obser- 
vational evidence from CR abundance ratios is -2.4. This value still results in a discrepancy 
from observation by a factor 10. 

3.2 Anisotropy and Galactic distribution of CRs 

Biermann and Davis also directed their attention to air-shower experiments that had just 
shown that the amplitude of anisotropy of Galactic cosmic rays was surprisingly small, 
below 10 -3 up to 10 eV and still below 10 % at 10 17 eV. These experimental results have 
stood the test of time[||, [94], |l], 

If cosmic rays are produced mainly within the solar circle (as is the case if SNRs are their 
dominant sources) at earth their diffusive leakage out of the Galaxy is equivalent to a slow 
net motion towards the Galactic anti center. A finite anisotropy is then expected due to 
the so called "Compton-Getting" effect fl32fl : less particles are incident in the direction of 



net motion. The most recent detailed analysis of this effect [86] finds that the amount of 
anisotropy expected at an energy of 10 14 eV is far larger than observed in three out of 
four considered models (by factors between 5 and 40). The remaining model manages to 
reproduce an anisotropy of 10 -3 - but only at the price of implausible ad hoc assumptions, 
like a CR source spectrum of the form ~ R _2 4 /(l+(2/R) 2 ) 1 / 2 (unexpected in any theory) 
- here R is the rigidity in GV - and an energy dependence of a cosmic-ray diffusion constant 
in flat contradiction to the most recent analyses of CR abundance ratios [|70||. 



13 This is only true of the CR producing SNRs are already known. The idea of an hitherto unknown class 
of "SNRs" with very different properties is indeed plausible and the subject of some proposals presented in 
section ^| 

14 Because the unity of the spectrum up to the knee implies a common origin of all CRs below the knee. 
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This problem is sometimes dismissed as mere coincidence. Perhaps the effects of various 
local sources of CRs conspire to result in an unusually small anisotropy near earth? Such a 
"conspiracy" appears unlikely when considering the determination of the CR density within 
the Galaxy by using observations from the 1970s onwards of low-energy 7-rays. The small 
anisotropy turned out to be the first indication of a deep lying problem. 
The radial scale length of the cosmic-ray density - derived in the interval between 5 and 
20 kpc - has a value between 16 kpc and 34 kpc in various recent determinations based 
in data obtained with the EGRET experiment |99|, |46| , ^(earlier determinations based on 
SAS and COS-B data had given similar results). The space-density of SNRs is maximal at 
a distance about 4 kpc from the Galactic center and falls off with a scale length of about 5 
kpc with increasing Galactocentric distance. Calculation of CR propagation based on this 
distribution predict a Galactocentric gradient which is about 3-5 times smaller than the 
observed one (see fig.||) . The distribution of CRs within the Galaxy is much more homoge- 
neous than expected, and this leads to the smaller than expected isotropy. There are two 
possible explanations for the discrepamcy. Either the Galactic-propagation properties or 
the source location of CRs are different from what was thought plausible up to now. 
With remarkable foresight Biermann and Davis had anticipated the problem just discussed 
from the limited data available to them and proposed an age of the observed cosmic rays of 
several billion years - one to two orders of magnitude larger than normally assumed. The 
longer diffusion time would indeed quantitatively result in the observed more homogenous 
distribution of CRs in the Galaxy. However, their proposal was ruled out by later exper- 
imental results: In particular, the work on the relatively large concentration of unstable 
isotopes with lifetimes in the order of 10 million years (Be 10 and Al 26 ) - performed after 1958 
- has convincingly proven that the CR ages cannot exceed several hundred million years. 
An alternative modification of the propagation properties suggests that the timescale of 
CR escape from the Galaxy might decrease towards the Galactic center pi]. This reduces 
the expected Galactocentric CR-density gradient. However the inferred total CR power of 
the Galaxy then exceeds the conventional value, thus requiring unreasonable acceleration 
efficiencies in SNRs (section |3.1|) . 

It is then most natural to blame the Galactic location of the CR sources for the small 
Galactocentric CR-density gradient. The main cosmic-ray sources are distributed 
in the Galaxy with a Galactocentric scale length in the order of 25 ± 10 kpc. 

Various systematic uncertainties are still too large to derive the specific form of the source 
distribution - it might e.g. be exponential with a scale length of order 25 kpc or, alterna- 
tively, constant up to 25 kpc with an additional Galactic-center contribution of order of 
10% of the total. 

The distribution of SNRs in the Galaxy has not only be observed directly (2^], but can also 
reliably be inferred - yielding similar results - from the one of their precursors, massive 
stars in our Galaxy and their brethren, pulsars^6|. This firmly excludes that SNRs may 
be distributed in the Galaxy as required by the assumption that they are the main sources 
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of CRs (as had been suggested by Strong and Moskalenko[101]). There simply is no known 



class of potential CR accelerators located as required. Even Galactic dark matter (if it 
exists) is too strongly concentrated towards the Galactic center. The main sources of 
hadronic cosmic rays must be a completely novel class astrophysical objects. 

3.3 The Origin of CR electrons in SNRs 

Biermann and Davis continue after their rejection of the hypothesis of SNRs as the main 
CR sources: 

"Supernovae might well maintain the energy and number balance of the relativistic electrons 
responsible for the radio radiation of the halos. This is independent of the overall energy 
balance of the nuclei, which constitute > 99 % of the cosmic rays. " 

I think that Biermann and Davis were correc t0- In SNR 1006 about 1% of the remnant's 
total kinetic energy was converted to relativistic electrons |69[][[]. If this value is typical for 
Galactic SNRs, diffusive and radiative losses of Galactic CR electrons can be replenished 
by SNRs. Electrons are observed to be accelerated in SNRs up to the highest energies 



observed in the ambient electron spectrum [103]. Moreover electrons - unlike the hadrons 



are not found to have a homogeneous distribution in the Galaxy. Rather they are strongly 
concentrated towards the Galactic center. The radial scale length of the synchrotron emis- 
sivity of the Galactic thick disk is only 3.9 ± 0.2 kpc - in good agreement with theoretical 
expectation under the assumption of an origin of the CR electrons in SNRs|J][5 
The excellent agreement between theoretical expectation and observational results under 
the assumption of a SNR origin of CR electrons makes the glaring failures of theory in the 
case of nucleons seem all the more alarming. 

4 Origin of "high-energy" (energies above the "knee") CRs 

4.1 Morrison's dark question: absence of a cutoff in the CR spectrum 
near the knee 



Philip Morrison|76| already realized that SNRs cannot accelerate CRs to energies above 
about 10 14 eV (the "knee" feature in the primary CR spectrum lies at a rather higher 
energy of 3 x 10 15 eV). His conclusion has since then been confirmed in numerous detailed 
theoretical studies |62[. Impressive observational evidence on the reality of the cutoff pll] 



15 "Novel" is meant in the sense "not yet directly observed" , the new source-class might well be related to 
known phenomena, like supernovae. 

16 See Dar et al.]35j for a different view. 
This value was derived from the expressions in Ref. |j9j and the experimental flux at TeV energies) 103] 



17 

18 r 



This conclusion only holds under the observationally plausible assumption that the mean Galactic 
interstellar B-field falls with a Galactocentric radial scale length much in excess 5 kpc. 
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PI make it unlikely that recent ingenious suggestions (e.g. fl2| ) to circumvent Morrison's 
limit are realized in SNRs. Morrison called the question of the origin of CRs with higher 
energies "dark", a designation that remains completely appropriate to the present dayp*|]. 

4.2 Today's second dark question: absence of the Greisen cutoff 

The most recent results from the AGASA air shower experiment indicate that the CR spec- 
trum continues with a differential index of about -2.5 from an "ankle" at 5 x 10 18 eV to 3 
x 10 20 eV without any evidence for a spectral cutoff above 5 x 10 19 eV that is expected if 
UHE CRs are extragalactic[102]. The newest data|93|] favor a power law from 5 x 10 18 eV 



to 3 x 10 20 eV without any structure. Three basic explanations have been brought forward 
to understand this(7|Q: 

1. Magnetic fields with a strength commonly encountered only in the center of large clus- 
ters unexpectedly exist also in extra-galactic space within about 10 Mpc of the Milky Way. 
This allows to "wash out" both anisotropics and the Greisen cutoff to a smooth increase 
of the spectral index in the mentioned energy range . UHE CRs could then be due to 
active galaxies @ 0. 

2. New physics occurs [|97|]. The ideas proposed have a similar status as the ones of Millikan 
and Lemaitre mentioned in the introduction. 



3. UHE CRs have their origin in the Galactic halo[^, 59, pq| . 

I concluded in section [T^ that the sources of low-energy CRs are a novel type of astro- 
physical object that is located mainly in the outer reaches of our Galaxy. The simplest 
point of departure for further research is then to adopt option 3 and to postulate that 
low and UHE CRs - and then indeed CR of all energies - have a common ori- 
gin in the same novel source class. A "numerical coincidence" is in favor of such a 
paradigm. Magnetic confinement to our Galaxy is usually expected to set in at energies 
below a few times 10 18 eV. The "ankle" is an increase in the spectral index of the primary 
CR spectrum by about 0.5 at about this energy. If only one source class dominates CR 
production at all energies, this coincidence can easily be naturally understood as the onset 
of Kraichnanian|28| Galactic confinement [105, p0[ . 



19 Reynolds and Keohane conclude that in 14 SNRs they studied, only one (none) has a maximum electron 
energy above 100 (1000) TeV and that this limit should also apply to nuclei. 

20 I do not list the possibility of an origin in the Galactic disk that had been resurrected by somejl^]. It 
requires an iron dominated composition above 10 19 eV that was experimentally ruled out recently 

21 In addition, a density ratio of CR producers inside and outside the local supergalaxy higher than the 
one suggested by independent astrophysical evidence has to be assumed |l8). It is obvious that this proposal 
was born of distress. 
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5 Alternatives to SNRs as CR major sources 



To my knowledge there have been only two proposals about an origin of the bulk of cosmic 
rays proposing objects different from SNRs or pulsarsPI within the last decade. In addition 
there is one historic proposal for an origin of CR in the Galactic halo. 

5.1 Burbidge 



In 1955 Geoffrey Burbidge[23] proposed that CRs are accelerated up to the highest energy 
in the Galactic halo. Density irregularities with typical velocities of 200 km/sec were 
proposed as scattering centers. The ultimate energy source was imagined to be Galactic 
rotation tapped via turbulent hydromagnetic frictional effects. 

5.2 Dar and Plaga 



Gamma-Ray Bursts (GRBs) are non-thermal objects that are ill understood [48]. We are 
sure that sometimes GRBs occur in our own Galaxy, but where, how often and with what 
kind of remnants remains shrouded in mystery. It is then an obvious speculation to identify 



the main source of CRs with GRBs (section |3.2[). This has been first proposed by Arnon 



Dar[34] and worked out by Dar and Plaga pq]. In the latter paper GRBs are proposed 
as "universal" (dominating at all energies) CR sources. An origin of CRs in GRBs (at 
energies above the "knee") was first proposed by Milgrom and Usov[71], Vietri|106, |107j ] 



and Waxman |lll| .) More recently Plaga[^2| proposed that bipolar supernovae[10£ , 27] eject 
plasma clouds into the Galactic halo. These "plasmoids" are posited as the main sources 
of CRs. He based this scenario on the "cannonball model" of Dar and de Rujula|37|, pSj] . 

5.3 Dermer 



It might be that Dar and Plaga pq] fell into the same trap as Shklovsky[95] and Burbidge [p^] 
before: again the source of CRs is directly identified with a bright non-thermal object in 
the sky. Charles Dermer recently suggested |4(| that - in a addition to GRBs - a special type 
of stellar collapse event resulting in "fireball transients" (FT) occurs in the Galaxy. FTs 
are explosion events where > 50 % of the kinetic energy of the ejecta resides in relativistic 
(r > 2) outflows. These outflows are proposed as the main source of Galactic CRs. GRBs 
would then only be the exceptionally baryon poor "tip of the iceberg" of the underlying 
novel FT source class, that is up to now mainly seen in the "light" of hadronic CRs. In 
Dermer's scenario UHE CRs are produced on extragalactic GRBs, whereas Dar and Plaga 
propose Galactic GRBs also for this component. 



22 The known population of Galactic pulsars falls short by about an order of magnitude in power to 
replenish the loss of Galactic CRs, even if high acceleration efficiencies are invoked. 
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5.4 General constraints on alternatives 

I expect that the final answer will involve two ingredients: an origin mainly (but probably 
not exclusively) in the Galactic halo (taken from Burbidge and Dar/Plaga) and a process 
with a more direct and effective conversion of gravitational into kinetic energy than the one 
that shaped SNRs (taken from Dermer and Dar/Plaga). This leaves room for ideas quite 
different from the ones presented in this section, like an origin in a novel, primordial class 
of neutron or quark stars in the Galactic halo. 

6 Conclusion - present status and future progress. A per- 
sonal view 

6.1 Present status - my view 

Summarizing, SNRs are not the main sources of Galactic cosmic rays. SNRs are the main 
sources of one of its minor components, the electrons. The distribution of the actual main 
sources of Galactic cosmic rays differs from that of any other known class of potential CR 
accelerators. It stands to reason to speculate that these sources accelerate the complete 
CR spectrum up to the highest energies. The astroparticle physics community holds in 
its hands solid evidence for a completely novel type of astrophysical phenomenon. It must 
rise to the challenge to study this phenomenon, not by merely waiting for a serendipitous 
observational discovery but by systematic theoretical work on alternatives to a SNR origin 
of CRs. 

6.2 Present status - prevailing view 

The extremely disappointing experimental situation is not denied in the astroparticle/cosmic- 
ray community. It is generally and clearly stated that a CR origin mainly in SNRs is not 
certain^]. However, u faute de mieutf 1 ("for lack of something better", quote from a report 
by a working group on "tests of Galactic cosmic-ray source models") any consideration 
of alternatives to SNRs as major CR sources is refused (not a single alternative to a SNR 
origin of the low-energy component is mentioned ^). Thus the prevailing view is that 
theoretical work on the origin of the low-energy component must remain restricted to the 
presently known class of SNRsp^l, because it is "the (only) one model that has been 
worked out in some detail\$3\"f% 

23 E.g. Matthew Baring skilfully avoids a personal endorsement in his summary of the 1999 26th ICRCj^] 
("... has led to the almost universal acclaim that SNRs are the site of CR acceleration. Yet ...") 

24 The strongly related concept of super bubbles - agglomerates of several old SNRs - as a possible site of 
CR origin ji?]] is discussed, but the authors stop short of considering them as the primary CR source class. 

25 formed of non-relativistic, very roughly spherical, SN ejecta 

26 In support of my evaluation I quote from the rejection of an anonymous referee on a manuscript 
suggesting an alternative to a SNR origin of CRs[p2j. The referee had demanded an impossibly large 
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6.3 If the prevailing view continues to predominate... 

...this will have two unhappy consequences: 

1. The problem of the origin of CRs will probably be solved outside the 
astroparticle-physics community, when eventually other astrophysicists will stum- 
ble on the CR source class by chance^]. 

2. The vital chance will be lost to be the first to study the new source class by 
means of CR evidence. 

To illustrate the latter point, Arnon Dar and I proposed an extremely jetted geometry of 
GRBs (later worked out in detail by him and Alvaro de Rujulaj3?J within the "cannonball 
model") guided by CR evidence [j36|, 

6.4 Future progress - some important open questions 

Two questions are of primary importance for the problem of CR origin and serve as a 



possible point of departure for the research recommended in section 6.1. The third one is 
of more general interest. 

1. There is solid observational evidence that some fraction of supernovae involves the 
ejection of ultrarelativistic plasma with a total kinetic energy larger than the one of 



SNRs |109| , [27fl (the most recent observation of a gamma-ray burst [ 48 1 - supernova asso- 
ciation is discussed by Lazzati et al. |]65| ). 

Where are their remnants (observation) and what properties do they have (the- 
ory) 0? What is the CR total -power and spectrum produced in these events 
in our and other galaxies? 



2. There is solid observational evidence |110| , p4j] that spiral galaxies possess the extended, 



thermal, magnetized halo or "corona" originally introduced by Spitzer|95], that our Galaxy 



forms no exception [113] and that these halos are in a very dynamic, turbulent state[110 



Where and with what total power are particles accelerated in these halos? 
3. What is the total cosmic-ray power of our Galaxy? 

If the answer turns out to be significantly higher than the standard value of 1.5 x 10 41 
erg/sec[^7| (as recently advocated by Dar & De Rujula[pS|]) a similar upscale is necessary 
for other galaxies. With an increase by more than a factor 10 CR pressure is expected to 
dominate the thermal gas pressure in the center of galaxy clusters [Bq]. This would refute 



extension of an already relatively long manuscript and reasoned: "Dr. Plaga argues that 15 A & A pages is 
quite a lot, but he is trying to challenge the cosmic ray edifice that stands upon uncounted thousands of pages, 
articles, and conferences F' . Moreover in none of two recent reviews of CR origin^, any alternative in the 
sense of section ^ is mentioned, though in both cases such alternatives were presented at the corresponding 
conferences. 

27 Because VHE astrophysics lacks sensitive all-sky detectors, it seems likely that the objects will even- 
tually be discovered in the X-ray. 

28 A consideration of various anomalous non-thermal remnants in our Galaxy |l00| , |Io| , |8iJ might be 
interesting in this connection. 
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one of the strongest arguments Q against "MOND", Mordehai Milgrom's modified Newto- 
nian dynamics [72]. 
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Figure 1: The flux of hadronic CRs rays as measured by various experiments as a function 
of total energy. The faint dashed line through all the data is a power law with a constant 
index of -3, and demonstrates that a single power-law is an excellent first-order approxi- 
mation to the experimental data over the whole energy range. The deviation from a power 
law at energies below 10 GeV is due to solar and terrestrial modulation effects. The basic 
components of the experimental CR spectrum in the "eclectic (many source classes)" sce- 
nario, generally accepted presently, are sketched. More than four decades have passed since 
P. Morrison suggested a slightly different "hierarchical" framework [76], in which cosmic- 
rays fill increasing shares of the universe with increasing energy. From the lowest energies 
up to the "knee" at 4 x 10 15 eV SNRs are the main contributors (full line). Theoretical 
arguments and experimental evidence suggest an earlier cutoff around 10 14 eV, this is sym- 
bolized together with the question marks. At higher energies up the the "ankle" at about 
5 x 10 18 eV another source class of unclear nature, labeled here with "Unknown sources" 
(dashed line) takes over. At still higher energies an extragalactic component (dash-dotted 
line) dominates, that is expected to cut off at the so-called "Greisen limit" of 5 x 10 19 
eV (the fact that the observed spectrum shows no evidence for such a cut-off is symbol- 
ized by the question marks). The missing cut-offs and exceptionally smooth joints of the 
three components (in particular at the "knee") finds no natural explanation at the moment. 
The idea that this theoretical patchwork does not do justice to the impressing simplicity 
and unity of the experimental data strongly suggests itself. It seems more reasonable to 
posit that the "Unknown sources" component dominates the hadronic CR luminosity at all 
energies. (Figure modified from an original of S.P.Swordy) 
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Figure 2: This X-ray image of the Cassiopeia A (Cas A) supernova remnant is the official 
first light image of the Chandra X-ray Observatory fjOf , The 5,000 second image was made 
with the Advanced CCD Imaging Spectrometer (ACIS). Two shock waves are visible: a 
fast outer shock and a slower inner shock. The inner shock wave is believed to be due 
to the collision of the ejecta from the supernova explosion with a circumstellar shell of 
material, heating it to a temperature of ten million degrees Celsius. The outer shock 
wave is progressing into the normal interstellar medium. The bright object near the center 
may be the long sought neutron star or black hole that remained after the explosion that 
produced Cas A. 
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Figure 3: The region around DR4 in the light of 92 cm radio waves (left image, from the 
WENSS survey |115fl ) and 100 //m infrared light from the IRAS ISSA server |6l[| . The radio 
image indicates the distribution of relativistic electrons emitting synchrotron radiation in 
the magnetic field of a molecular cloud compressed by the SNR shock wave. The infrared 
emission traces dust in the cloud heated by the passage of the shock wave. The similar 
distribution of both is an argument in Cavour of an interaction of shock wave and cloud. 
Both frames are untreated 1-degree cutouts from the respective cataloguers. They are 
centered on the coordinates ra: 20h 22m 37.98s; deck +40d 09m 41.0s (J 2000). 



25 




Figure 4: Various experimental upper limits on the flux of VHE 7-rays from the region 
DR4 in the SNR G78. 2+2.1 are compared to a theoretical prediction, derived under the 
assumption that SNRs are the main sources of hadronic rays (from the thesis of Proschjp^] 
which contains all the references.) The predictions under the assumption of two spectral 
power-law indices of the accelerated hadrons -2.1 (full line) and -2.3 (dashed line) are shown. 
The universal value for strong shocks theoretically expected is -2.1. A spectral cutoff at 
100 TeV was also assumed. The detailed assumptions used for deriving the theoretical 



curves are explained in Prosch et al.[84|. This publication contributed the point labeled 
"AIROBICC 1996" and was the first to officially announce the failure of the theoretical 



predictions by Aharonian et al.[42|. 
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Galac to centric radius, kpc 



Figure 5: Radial distribution of 3 GeV protons at z = 0, for diffusive re acceleration model 
with different halo sizes Zh = 1, 3, 5, 10, 15, and 20 kpc (five continuous solid curves). The 
source distribution assumed for the calculation of these curves is that for SNRs as given by 
Case & Bhattacharya 25], shown as a dashed line. The cosmic-ray distribution deduced 
from EGRET >100 MeV gamma rays (|99[|) is shown as the histogram. It is seen that the 
distribution of hadronic cosmic rays is essentially homogenous up to about 15 kpc. If cosmic 
rays are mainly produced in SNRs or similarly distributed objects a strong concentration 
towards the Galactic center is expected - in strong contradiction to the data. (Figure taken 
from Strong and Moskalenko II Oil].) 
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